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The form factor and the couphng constant in the r]i,BB* vertex are evaluated in 
the framework of three-point QCD sum rules. The correlation functions responsible 
for the vertex is evaluated considering contributions of both B and B* mesons as 
CN . off-shell states. The form factors obtained are different for the two cases, whereas 

■ the final results of the coupling constant are compatible. 
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I. INTRODUCTION 

the bottomonium ground state rji, was discovered in the processes T{3S, 2S) — t- 
^ . , , |2[. More data are expected to be accumulated in the forthcoming SuperB and 

QJ . LHC-b. Therefore, investigation of the spectroscopy and decay processes related to this 

state becomes rather instructive. Since it is still difficult to study strong interaction phe- 
nomena at non-perturbative regime, the study of quarkonium decay is generally performed 
^ ■ in the framework of effective Lagrangian with meson exchange. Thus, the reliable determi- 

50 ■ nation of various characteristics, such as form factors and coupling constants are needed. In 

^ ■ the decay process T(45', 5S) -> ?7b + 7, rescattering effects are important for understanding 

O ■ the anomalous largeness of the branching ratios [3| . As suggested in their work, this decay 

^ ; proceeds in two steps. First the T(4S', 55) decays into a BB* intermediate state, and then 

these two particles produce the final states rjb and 7 by exchanging a B meson. In order 
to compute the effect of these interactions in the final decay rate, the coupling constant 
^ i of the r]j,BB* vertex is an necessary input parameter. To describe strong interactions of 

r]bBB* at the hadronic level, the following effective Lagrangian is employed [sl: 



^ViBB* = ig^,BB'B*d''f]bB^. (1) 

The rjbBB* interactions are characterized by strong coupling constant g-q^BB*- However, 
such low-energy hadron interaction lie in a region which is very far away from the per- 
turbative regime, precluding us to use the perturbative approach with the fundamental 
QCD Lagrangian. Therefore, we need some non-perturbative approaches, such as QCD 
sum rules [4|-[6j, to calculate the form factors. 

In this article, we calculate the form factor gr)tBB*{Q'^) in the framework of three-point 
QCD sum rules (QCDSR). More specifically, we evaluate two form factors: one when B is 
the off-shell particle and another when B* is off-shell. The two results are parametrized by 
analytical forms and then extrapolated to obtain the coupling constant g-q^^BB*- Herein, we 
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use the same technique for the study of the couphngs in the vertices D*D7i [7|, l8|, DDp 



D*D 



D*D*p DDuj 



PJJ, D*DK*{892) [13], DsDK* p4 and B*^B*K (15|. 
The outhne of the Letter is as follows. In Sec. [Tll we present QCD sum rules for the 
considered vertex when both B and B* mesons are off-shell. Sec. IIIII is devoted to the 
numerical analysis and discussion. 



II. THE SUM RULE FOR THE r]bBB* VERTEX 

In this section, we present QCD sum rules for the form factor of the rji^BB* vertex. The 
three-point function associated with the rj^BB* vertex, for an off-shell B meson, is given 
by 

rj(p,p') = / d'xdS e'^'-e-'^-y{0\T{j^\x)f^{y)f^^mO), (2) 

where the interpolating currents are j^* {x) = q{x)'y^b{x), j^{x) = iq{x)'~^^h{x) , and 
j'^''{x) = ib{x)'-f5b{x) . The correlation function for an off-shell B* meson is 

rf = / d'xd'y (^P'-e-^^-y {0\T{fix)jP{y)f''^iOm, (3) 



and q = p' — p is transferred momentum. The general expression for the vertices ([2]) and 
(IH]) can be written in terms of the invariant amplitudes associated with two independent 
Lorentz structures: 

T,ip,p') = T,ip',p'\q')p, + Up^p'\q')p'^. (4) 

In compliance with the QCDSR, the above correlation functions need to be calculated 
in two different ways: in the theoretical side, they are evaluated by the help of the operator 
product expansion (OPE), where the short and large distance effects are separated; In the 
phenomenological side, they are calculated in terms of hadronic parameters such as masses, 
leptonic decay constants and form factors. The sum rules for the form factors are obtained 
with both representations being matched, invoking the quark-hadron duality and equating 
the coefficient of a sufficient structure from both sides of the correlation functions. To 
improve the matching between the two representations, double Borel transformation with 
respect to the variables, = —p^ — )■ and P'^ = — p'^ — ?■ M'^, is performed. 

The physical side of correlation function ([2]) is studied with hadronic degrees of freedom. 
From Eq. ([T]), one can deduce the matrix elements associated with the rjbBB* vertex: 

{B{p)7^b{q)\B*{p + q,e)) = g,,BB* {p ■ e*) ■ (5) 

The meson decay constants /^j, Jb, and Jb* are defined by the following matrix elements: 

2 



{Q\f\B{p)) = ^fB, 
nib 
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(0|j« \B%p)) = TUB^fB^eUp)- (6) 
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FIG. 1: (a) and (b): Bare loop diagrams for the B and B* off-shell, respectively 



Saturating Eq. ([2]) by the complete set of appropriate B*, B* and r^b states, using Eqs. ([5 
and on]), and then summing over polarization vectors via 



2 ' 



the physical side of the correlation function for B off-shell is obtained 

C 



+. 



In a similar way, the physical side for an off-shell B* meson is obtained as 

C 



X 



B) 



9iBB*{<l'){-9^u + ^)p, 



m 



B* 



+. 



where we use the abbreviation 



C 



mlm\mB*f^^jB*fB 



2ml 



(7) 



(8) 



(9) 



and represents the contribution of the higher states and continuum. 

In the following, we turn to the computation of the correlation functions in the QCD 
side. Each invariant amplitude r*(p',p) appearing in Eqs. (jl]) can be written in terms of 
perturbative and condensate terms 



(5) 



(10) 



where T^^^ is the perturbative contribution, and ■ ■ -, vf' are contributions of conden- 
sates of dimension 3, 4, 5, • • • operators in the OPE. The perturbative contribution and 
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FIG. 2: Diagrams for contributions of bi-gluon operator in the case B off-shell. 
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FIG. 3: (a) and (b): Diagrams corresponding to quark condensate for the B off-shell and B* 
off-shell, respectively. 



gluon condensate contribution can be defined in terms of double dispersion integral as 

1 
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where p^'^^{s,u,Q'^) is the perturbative spectral density and t = q^. The perturbative 
spectral density can be obtained by calculating diagrams (a) and (b) in Fig.([l]). In 
the calculation, Cutkosky rules are adopted to deal with the usual Feynman integral 
of these diagrams, i.e., by replacing the quark propagators with Dirac delta function 
2 (— 27ri)(5(g^ — rn?)6{q^). The physical region in s and u plane is described by the 



following inequalities: 
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where A(a, b, c) = + 6^ + — 2ac — 2bc — 2ab and t = = —Q^. The diagrams for 
the contribution of the gluon condensate in the case B off-shell are depicted in Fig. ([2]). 



We follow the method employed in Refs. la, Il7|], namely, directly calculate the imaginary 



part of the integrals in terms of the Cutkoskys rule. In this Letter we use the structures 
for the off-shell B meson and p'^ for the off-shell B* meson. After some straightforward 
calculations, the spectral densities is obtained as following: 

pf^^'"'^(s,t,M) = — - — ^ ,^i3/o (3 * {u{2t{ml - mbrriq - s) + s{mb - mqf) 
+ (mb - mq){s - t){mi,t + nigs) + mbU^{mq - mi,))), 
= wSt^M^H-*- 5^). (12) 

for the off-shell B meson, and 

^^'''■^(s, t, u) = J- — ^ * {ml{-{s - t + u)) + ml{s - t + u) + u{s + t - u)), 

for the off-shell B* meson. All powers of the light quark mass are included when calculating 
the spectral density. In the above expressions, pi(2) represents the spectral density for 

n □ 

As what has been shown in p, 17|, heavy quark condensate contributions are negligible 
in comparison with the perturbative one. Thus, only light quark condensate contributions 
to the correlators are calculated with considering the diagrams (a) and (b) of Fig. ([3]). 
However, contributions of the light quark condensate contributions are zero after the double 
Borel transformation with respect to the both variables and P'^. The D = 5 quark- 
gluon mixing condensate contributions are also zero. 

Quark-hadron duality assumption is adopted to subtract the contributions of the higher 
states and continuum, i.e., it is assumed that 

phigherstates^^^ m) = p^^^{s, U, t)e{s - So)e{u - Uq) , (14) 

where Sq and Uq are the continuum thresholds. 

The double Borel transformation with respect to P^ = —p^ — )■ and P'^ = —p''^ — )■ 

M'^ are applied when matching two sides of the correlation function. The final sum rules 
for the corresponding form factors are obtained as: 



5(1) .n2^ (Q^ + ^ 4^ 
9r„BB*{Q ) = ^ e^e^ 



so ruo 

ds / du 



4 Vr^ J 4ml J ml 

h b 

(pf ^''■^ (s, t, u) + pf (s, t, u))e[l - M)2]e^e^ 



(15) 
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and 



so ruo 

ds I du 



4 TT^ Jiml Jml 

b b 



2771%, + ml.* ) '"'76 Hb 
{ml - ml -t)C 

{pT^''''\s, t, u) + pT^'\s, t, uMl - F^' {s, ufje^e^F 



(16) 



It is noted that in the following analysis, we use the relations between the Borel masses 



^2 , .„ ^2 



M2 and M'^ as = ^ for a S off-shell and ^ = ^ for a B* off-shell. 



III. NUMERICAL ANALYSIS 



In the numerical analysis of the sum rules, input parameters are shown in Table [11 We 
take and /^^ from Ref. The continuum thresholds, sq and mq, are not completely 
arbitrary as they are correlated to the energy of the first excited states with the same 
quantum numbers as the states we concern. They are given by Sq = {iTirn, + ^s)^ and 
Uq = {m + 1 where m is the B* meson mass for the case that B is off-shell and the B 
meson mass for that B* is off-shell. A„ and are usually around 0.5 GeV. The threshold 
So, Mo and Borel variable are varied to find the optimal stability window where OPE 
convergence, stability and pole dominance of the sum rule with the Borel mass parameter 
are satisfied. 

TABLE I: Parameters used in the calculation. 

m,(GeV) mi,(GeV) mi,,(GeV) m^,(GeV) /^(GeV) [19j /^^(GeV) [20] /^.(GeV) 
4.2 5.28 5.325 9.4 0.21 0.16 0.705 



Using As = A„ = 0.5 GeV for the continuum thresholds and fixing = 2GeV^, the 
dependence of the S'^bb* (<5^ = 2GeV^) on Borel mass are shown in Fig. (jlK)- From this 
figure, we see the results exhibit a good OPE convergence and stability for > 20 GeV^. 
Fig. (jib) demonstrates the contribution from the pole term with variation of the Borel mass 
M^. We see that the pole contribution is larger than continuum one for < 35GeV^. 
We choose = 33 GeV^ as a reference point. 

Now, we would like to discuss the behavior of the form factors in terms of Q^, which 
are plotted in Fig. [5l In Fig. [5l the circles correspond to the form factor g^^BB*{Q'^) the 
interval where the sum rule is valid. For g^^BB*{Q'^)i result is better extrapolated by 
the Gaussian fit parametrization, 

rn2^ 9n9P ^ -(Q^ + 21.9GeV^)^ 

The coupling constant is defined as the value of the form factor at = — m^, where m 
is the mass of the off-shell meson. Using = —m^ in Eg. (1171) . the coupling constant is 
obtained as g^j^BS* 

= 19.0. 
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FIG. 4: a) The dependence of the form factor g^jj^^* 

(Q2 = 2.0 GeV^) on Borel mass for 
As = 0.5 GeV and A„ = 0.5 GeV. The notations a, /? and 7 correspond to total, perturbative 
and four-quark condensate contributions respectively and b) pole-continuum contributions. 
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FIG. 5: g^^BB-iQ"^) (circles) QCDSR form factors function of Q^. The solid lines correspond 
to the Gaussian parametrizations. 
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FIG. 6: a) The dependence of the form factor g^^BB'iQ'^ ~ 2.0 GeV^) on Borel mass parameters 
for As = 0.5 GeV and = 0.5 GeV. The notations a, /? and 7 correspond to total, perturba- 
tive and four-quark condensate contributions respectively and b) pole-continuum contributions. 



In the case that B* is off-shell, Fig. (jS^) demonstrates a good OPE convergence and 
stability of g^^BB* with respect to the variation of Borel mass for > 30GeV^. We 
see that the pole contribution is bigger than the continuum one when the Borel mass 
M2 < SOGeV^ from Fig. ([6b). Fixing = 40 GeV^ our numerical results can be fitted 
by the Gaussian fit parametrization 



^3*(Q') = 26.1 Exp 



-(g2 + 30.7 GeV' 
951.6 GeV^ 



2\2 



shown by the solid line in Fig.©. The coupling constant g^*BB* 



-m 



B* 



in Eq. 



25.8 is obtained at 
(fT5]l . The two parametrization processes lead to the compatible 



coupling constant, which is a check on the extrapolation procedure. Taking the average of 
the two results, we get 



g^,BB* = (22.4 ±3.4). 



(19) 



In order to study the dependence of our results with the continuum threshold, we vary 
As u in the range 0.4 GeV < Ag^^ < 0.60 GeV. As can be seen in Fig. |H1 this procedure 
gives us an uncertainty interval of 13.8 < gr)tBB* 

< 26.4 and 17.7 < g^^'j% < 34.2. Taking 
into account the uncertainties due to the continuum threshold parameter, we finally obtain 



g^^BB* = 24.0 ± 10.2. 



(20) 



We noticed that, in Refs. 



22 



the value gr^^BB* = gr{is)BB = 15 is used as an input 
parameter. Their estimate is based on the assumption of the heavy quark spin symmetry. 
Considering the uncertainties, our results are compatible with their estimates. However, 
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FIG. 7: g^^BB*iQ ) (circles) QCDSR form factors as a function of Q . The solid lines correspond 
to the Gaussian parametrizations. 

our calculation shows that the predicted central value from QCDSR is 60% larger than 
their result. This may be a demonstration of some degree of heavy quark spin symmetry 
violation. The numerical value of this parameter may affect the physical observable to some 
extent. Thus, we expect a further theoretical investigation may give a better confirmation. 

In conclusion, we have used three-point QCDSR to calculate the form factor of the 
rjbBB* vertex. Both cases that B is off-shell and B* is off-shell have been considered. As 
a side product of the form factor, the coupling constant is estimated. 
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